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ABSTRACT 


Considering relative misalignment between the seal faces to be the 
basis of pressure generation and sealing, it is investigated 
whether hydrodynamic and squeeze film pressures provide stable 
seal operation. The motion of the flexibly moimted stator ring 
in response to the rotation of the rotor ring consists of a 
wobble. It is described by two non-linear differential equations 
for two mutually perpendicular components of the net tilt of the 
stator. The equations are reduced to linear form by assuming 
constant wobble frequency and conditions are found for the tilt to 
increase, decrease or remain constant. A relationship between the 
wobble frequency of the stator ring and the frequency of rotation 
of the rotor ring in terms of the support stiffness, moment of 
inertia and mean radius of the stator ring is established. 

It is found that when the condition for constant tilt is 
satisfied the wobble frequency of the stator ring is one-half of 
the shaft or rotating ring frequency. This state results in sero 
fluid film moments, causing failure of the seal. The relationship 
between frequency of rotation of the rotor and the support 
stiffness, moment of inertia and mean radius of the seal provide a 
means for the designer to avoid this mode of seal failure. 
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Considering relative misalignment between the seal faces to be the 
basis of pressure generation and sealing, it is investigated 
whether hydrod7n^uaic and squeeze film pressures provide stable 
seal operation. The motion of the flexibly mounted stator ring 
in response to the rotation of the rotor ring consists of a 
wobble. It is described by two non-linear differential eqixations 
for two mutually perpendicular components of the net tilt of the 
stator. The equations are reduced to linear form by assxaaing 
constant wobble frequency and conditions are found for the tilt to 
increase, decrease or remain constant. A relationship between the 
wobble frequency of the stator ring and the frequency of rotation 
of the rotor ring in terms of the support stiffness, moment of 
inertia and mean radius of the stator ring is established. 

It is found that when the condition for constant tilt is 
satisfied the wobble frequency of the stator ring is one -half of 
the shaft or rotating ring frequency. This state results in zero 
fluid film moments, causing failure of the seal. The relationship 
between frequency of rotation of the rotor and the support 
stiffness, moment of ineirtia and mean radius of the seal provide a 
means for the designer to avoid this mode of seal failure. 
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The critical condition to avoid is w = C 



m. where r_ 


is the mean radius of the seal ring, its transverse moment of 
inertia , k a measure of its supx>orb stiffness and c a constant of 


proportionality determined by the manner in which the support 


stiffness is modelled. 
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NOMSNCLATDRE 


h local film thickness; seal gap height 

K mean film thickness 


H 

I 

k 

1 

M 

0 

P 

r 

t 

X,Y,Z 

x,y,z 

r 

r 

A 

6 

e 

0 

P 


axial velocity due to change In mean film thickness 
moment of inertia 

unit vectors in the directions of the positive x.y and z 
axes . 

spring stiffness 
iincompressed length of spring 
moment 

geometric centre of the stator ring 
pressure 

radius; or cylindrical coordinate 
time 

inertial reference frame axes 
rotating reference frame axes 

relative angular misalignment (tilt) between seal faces 

tilting angular velocity 

spring displacement 

initial spring compression 

tilt parameter, yr_/H 

cylindrical coordinate; angular displacement , around seal 
ring measured from the positive y axis 


viscosity 



angular displacement around seal ring measured from the 
positive Y axis 

y/ Euler angle specifying the position of the x,y axes with 

reference to the X,Y axes in the horizontal plane 
r wobble freauency 

ci> rotor (shaft) angular velocity 

w absolute angular velocity of the x-y-z rotating 

reference frame 

5 absolute angular velocity of the stator ring 

Subscripts 

hyd hydrodynamic 

i inner radius 

m mean radius 

o outer radius 

\ 

3 sealed fluid 

sp spring 

X,Y,Z inertial reference frame axes 

x,y,z rotating reference frame axes 

xx,yy,zz indicate principal axes of the stator ring 
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1.1 Introduction 

A large number of seals of different types are used in 
Industry, for fluid sealing in systems with rotating shafts. The 
mechanical seal, as compared to simpler stuffing boxes and gland 
packings used for rotating shafts, has the advantage of smaller 
leakage losses, reduced maintenance and longer life due to no 
ware. The function of a seal is to separate pressurised fluids 
where, for instance, a moving shaft passes though a machine 
housing or passes from one part of the machine to another. It may 
also serve to prevent entry of foreign bodies into an operating 
medium or the loss of lubricant from bearings and transmissions. 
A typical low-cost volume-production face seal is illustrated in 
Figure 1.1. This type of seal is used in such applications as 
motorcar water pumps, domestic washing machines, and small 
centrifugal pumps. Very often the floating member is stationary, 
running against the end of the pump impeller for instance. The 
high pressure region is on the outside of the seal, causing a 
leakage path as indicated, "p^e mechanical seal restricts leakage 
of fluid in the radial direction, between the nominally plane 
sealing surfaces of the stator ring and the impeller. The 
floating stator ring is typically of a moulded resin material. The 
case is a one-piece pressed - steel ring within which a rubber 
bellows is fitted and retained by its own initial interference or 
a combination of this and the axial thrust of the face-loading 
spring. 
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rubber bellows 
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Fi£. 1.1 A low-cost volume-production mechanical seal 
(Crane Packin£ Ltd., Type 25A) 
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Mechanical seals, also known as radial face seals, have manv 
applications. In the engineering Industry, they serve to seal 
crank shafts, water pumps and auxiliary pumps in automobiles, 
power generating systems auad marine engines. They are also common 
in water turbines, boiler feed pumps compressors and high pressure 
rotary pumps. In aircraft and rockets, they are used for sealing 
gas turbine shafts, turbo-super chargers, hydraulic units and 
booster pumps for fuel and liquid gas. The fluids sealed include 
liauidoxygen, super heated water, solutions such as brewing 
liquors, dyes, sludge, various acids, hydro -carbons , asphalts, 
cocoa, plastics, radioactive gases and liquids. For various 
applications, suitable seal materials and special designs are used 
for chemical, mechanical and thermal compatibility. Pressures 
ranging fr<Mft high vacuum (lO”*' torr upto 200 bar) shaft speeds 
upto 60,000 rpm and temperatures from -200* C upto 450* C are some 
of the requirements. 

The stability of the seal operation is a crucial factor in 
order that the seal faces may not come into direct contact. Such 
contact can lead to frictional overheating and severe thermal 
distortion of the seal faces following which the seal falls. 

In addition, the following factors have an important 
influence on leakage losses, life and reliability. 

1. The width of the seal. 

2. The sliding speed. 
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3. The surface roughness and the waviness of the sliding 
surfaces . 

4. The temperature of the sealed liquid and the sliding surfaces, 

especially any variation with time. 

5. The shape of the interface profile. This is dependent upon 
mechanical and thermal distortions which can occur during 
operation. 

6. The choice of seal face materials. 

7. The sealed medium, its degree of solid contamination and its 
chemical behaviour. 

8. Oscillations, side loads, interrupted operation, periodic dry 
running, heating and cooling cycles, the direction of the 
radial flow leakage, eccentric rxinning and other factors 
determined by design on operating conditions, especially those 
affecting heat transfer characteristics. 

In Mechanical seals used under high loads, or where the 
sealing fluid has an inadequate lubricating ability, pressurized 
fluid is pumped into the seal interface. 

Successful operation of the mechanical seal in order to 
reduce wear and maintain the integrity of the sealing surfaces, 
primarily requires that the faces must remain separated by a 
lubricating film thoughout seal operation. Also, to minimize 
leakage, the seal - gap must be extremely small, of the order of 
10 mij . 
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Reference [13 shows that a lubricating film can exist between 
the seal faces. The paper also provides an insight into the 
mechanisms of successful seal operation. The mechanisms 
responsible for the development of the lubricating film pressure 
that acts to separate the faces of the seal, include angular 
misalignment, surface waviness including thermal deformation 
effects, surface asperities and axial vibration. 

Stanghan - Batch and Iny [2] applied the hydrodynamic theory 


of lubrication to 

calculate film pressure 

in 

an 

idealized 

wavy 

faced seal. The 

pressures were found 

to 

be 

sufficient 

for 


supporting the applied face load. On the basis of their results, 
Stanghan - Batch and Iny postulated that successful operation of a 
face seal depends on the generation of a hydrodynamic film 
between the seal faces. The film not only prevents direct contact 
between the two faces, but also provides the sealing pressure. 
The authors went on to suggest that the performance of mechanical 
seals can be improved by giving one face an appropriately wavy 
surface at the time of manufacture. 

An angular misalignment (tilt) is akin to a single lobed 
waviness . 

In this thesis, relative angular misalignment of the seal 
faces is considered as a possible mechanism for fluid film 
lubrication and for the development of the lubricating film 


5 




Fig. 1.2. Hydrodynaittic filn pressure genera'tion due "to relative 
face misalignment £2] 


6 



pressure at the seal interface. Equations for film pressures 
dependent on system parameters and seal geometry are derived. A 
model of the mechanical seal is considered, with one face rigidly 
mounted and the other with a floating and flexible movinting, with 
the object of analysing the dynamics and stability of the seal. 

Stanghan - Batch and Iny’s theory of face seal operation [2], 
envisages only a hydrodynamic film and hydrodynamic pressure. 
Since the hydrodynamic pressure is generated by a fluid film of 
varying thickness around the seal circumference, a 
hydrodynamically induced moment acts on the seal face. The 
effect of this moment is to neutralize the relative face 
inclination if such inclination consists of a single - lobed 
waviness of the seal face. A two - lobed waviness, on the other 
hand, generates a two-lobed hydrodynamic pressure which can 
produce a net zero moment and sustain the film conditions. Thus 
Stanghan - Batch and Iny were for c ed to suggest that the 
hydrodynamic action in a face seal can work only when a two - 
lobed waviness exists, ruling out the possibility of a relative 
angular misalignment as the origin of the hydrodynamic pressure. 

Figure 1.2 illustrates the action of the hydrodynamic moment 
and shows how it neutralizes the face tilt. In practice, the seal 
faces are carefully lapped before installation. The flexible 
movinting of one face of the seal however permits it to tilt, 
thereby generating a film with a varying thickness (single - 
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lobed wave). Any theory of face-seal operation should consider 
the possibility of hydrodynamic effects arising out of a relative 
angular misalignment and the conditions for the dynamic stability 
of such a configuration. 

The couple generated by pure hydrodynamic pressure tends to 
neutralize the relative tilt and establish parallelity of the two 
faces. In this process, fluid is squeezed out of the thicker 

regions of the fluid film. That, in turn, generates a squeeze 

film pressure which resists the hydrodynamic moment. Thus if 

hydrodynamic pressure can disturb the tilt configuration the 
squeeze film pressure may restore it. 

The objective of this thesis was to investigate whether 
hydrodynamically generated film pressure and squeeze film pressure 
provide stable seal operation. Therefore a simplified model of a 
mechanical seal is chosen, considering only the relative 
misalignment between the seal faces as the mechanism of pressure 
generation and sealing. A more complete model would include 

surface - waviness, axial vibration and the effect of the support 

arrangements. For axial vibration to take place the mean film 
thickness has to be time dependent. Support arrangement effects 
would have to take into account the flexibility of the seal faces 
due to the presence of rubber gaskets and 0-rings. The damping 
characteristics of the support arrangement would also have to be 
considered. In this first study, such effects are neglected. 
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Ft£. 1.3 Mathematical model of the face seal 
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1-2 The MechaPilcal Face Seal Model 

A simple mechanical face seal configuration in shown in figure 
1.1. Here the rotating seat is rigidly mounted on the shaft, with 
no axial translating freedom and with perfect alignment. 

The stationary ring (stator) is flexibly mounted on a spring. 
Therefore the stator can undergo both axial translation and tilt. 
The tilt can be looked upon as two rotations about mutually 
perpendicular diameters of the stator. 

Both faces are flat. The shaft has constant angular velocity. 
The lubricating film pressure that acts to separate the faces of 
the seal is developed because of the seal ring angular 
misalignment. The mathematical model of the mechanical face seal 
with a misaligned seal ring is shown in Figure 1.3. The seal ring 
angular misalignment creates a convergent and divergent region in 
the intervening fluid film. 

The principle of hydrodynamic lubrication is responsible for 
the film pressure development. 

In addition to the tangential motion causing hydrodynamic 
lubrication, the stator - ring can also have a motion normal to 
the surface of the rotor. This motion also helps in development 
of a film pressure and the principle is known as s<iueese film 
lubrication . 
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1 . 3 The Development of Fll« Presgmres 


The variation of film pressure in a mechanical seal is described 
by the Reynolds equation. The following assumptions are inherent 
in the use of the Reynolds equation: 

1) the fluid is incompressible and the flow is laminar, 

2) the fluid properties remain constant i.e. effects due to 
variation in temperature and pressure are neglected, 

3) inertial effects are negligible, 

4) the solid bodies remain rigid, 

5) the film is of sufficiently small thickness that the fluid 
pressure can be considered constant though the thickness of 
the film, 

6) the radius ratio r,/r (fig. 1.4) is of the order 1, i.e. the 

1 o 

face width is small compared to the mean radius of the seal. 
The simplified Reynolds equation for short beaurings, using 
cylindrical co-ordinates is, 

I?’ = ® 

The film thickness h is given by (refer Fig. 1.4), 

h = H + rr cos e (1.2) 
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where ^ is measured from the location of the maximum film 


thickness . 


Then, 


dh _ 

w - 


~rv sin e 


(1.3) 


If the film thickness is measured in a moving reference frame, 
then & is zero, as the time derivative of theta t.e. & is the 
angular velocity of the reference frame which will not be seen in 
the same reference frame. 

Therefore , 


^ = h -»■ yr cos O (1.4) 

Substituting equations (1,3) and (1.4) in (1.1), 

^ ~§f} ~ ~6^i«rr sin© + lE/^h + 12iJ rr cos © 

i.e 

r (“fe © « + 12p (H+ frcos ©) 

(1.5) 

On the right hand side of equation (1.5), the first term 
represents hydrodynamic effects due to « and the second tewn 

represents squeeze effect due to r and H. 


13 



1.4 Hydrodynamic Ef f<Kit 


As equation (1.5) is linear, the pressure distribution for only 
the hydrodjmamic effect is given by 

r ( 7 ^ “ -Bw'x sin O « 

“§? If) " " sin ^ 

Integrating once with respect to *r’ 

_ a 

(h*r ^ = -6^^x sin 6 — 

2 C 

4 dp _ -2fjcir sin © r a 

^•® dr ■ V® 1 .® 

h hr 

Substituting for h from (1.2) into ( 1 . 8 ) 

2 C 

dp _ sin dr. * 

_ 

(h + yr cos d) (h + yr cos d) r 


( 1 . 6 ) 

(1.7) 


2 C' 

_ - 2 M<«y sin d r ^ 

- 5. (i^aiooser,, 

E 

Integrating with respect to 'r’ using tables to integration [33 


r(i + r.SP.ggI)’ 
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-2M<«y stv0 h 

^ ~ tr» ■ •a 

n r cos ® 


(1 + »:5S2^) + 2 


^ ^rcondx - g(^^ >-coa9r ). 


i ^ 2-t- rcoser/H ^z ^ (____L__) 

* 1+ rcos©r/i5 1+ rcos^r/E 


+ C 


where C, end C, ere oonstents of lnte«retlon 
Simplifyinfi, 


^ _ -2tJ03r slug 

^ “ 9 9^ 

r cos & 


In (1+ 


yrcQsg ' ^ 




yvcos© * 


(l-2rr cose/h) 




yco3^/2h ^2 
1+ rcos0r/h. 


+ In (- 


■) 


1+ rooa0x/K 


+ c_ 


(1.9) 


In order to find the constants of tntesratton and C^, boundary 
conditions have to be utilized, at the inner and outer radius of 
the seal. 

The choice of having the sealed fluid placed on the inside of 
the sealing ring for an ‘inside-seal* or on the outside of the 
sealing ring for an 'outside-seal* , should be made to suit both 
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seal-leak&se, as well as duty reauirements . The disposition of 
the sealed fluid on the inside or outside of the seal is likely to 
have a strong impact on seal performance. Centrifugal action due 
to seal rotation tends to fling the interfacial fluid outwards. 
The radial variation of the hydrostatic pressure tends to drive 
the fluid along that pressure gradient, outwards or inwards as the 
case may be. Over and above these effects, there is the often 
observed phenomenon of "inwards pumping" that occurs in face 
seals. While inward pumping remains unexplained, Stanghan - Batch 
and Iny [2] suggest that the curvature of the seal faces 
produces a radial asymmetry in the pressure wave and induces a 
net inward flow of fluid across the faces. If the sealed fluid is 
placed on the inside of the sealing ring, this inward pumping 
effect can be utilized to suppress the flow of fluid acjross the 
faces due to the hydrostatic pressure gradient arising out of the 
centrifugal action of the sealed fluid pressure. However, this 
precise balance of flow is difficult to maintain in practice and 
should conditions change so that the inward pumping action exceeds 
the sealed fluid pressure, all fluid would be pumped out of the 
clearance between the faces, with disastrous consewences . Seals 
should therefore be designed to leak slightly to ensure that a 
lubricating film is maintained between the faces. The choice of 
an inside or outside seal configuration should also be made to 
suit the duty required of the seal. For moderate to high sealed 
fluid pressure, an inside seal should normally be used so that the 
inward pumping action can be utilized to limit the leakage. For 
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Inside seals with low sealed fluid pressures, and particularly 
with high rotational speeds, the inward pumping action may be two 
great to be overcome by the hydrostatic fluid pressure at 
practical clearances. A more satisfactory design would result 
from arranging the sealed fluid on the outside of the seal ring, 
in which case a steady supply of fluid to the seal interface is 
assured although the leakage rate may be high. 

In the present analysis an inside seal is chosen and the 
boundary conditions are p = p at r = r and p = 0 at r = r . As 

W ^ W 

equation (1.8) is much to cumbersome, further development of the 
pressure equation is postponed till after the develoiment of an 
approximated equation for hydrodynamic pressure. 

1.5 Sqaeeae Effect 

From equation (1.5), the pressure distribution for only the 
squeeze effect is given by 

“I" “li If) = 12 P (h + fr cosd) (1.10) 

ie §f) = 12 fi hr + 12 pyr* cose 

Integrating with respect to ‘r’ 

h r ” 12Mh — + 12My cosS — ^ 
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te 


_ * 2 C 

^ y 4. !L 

^ h* h* h®r 


Substituting for h from equation ( 2 ) 


dp 

dr 


QuEr 


(h+rr COS0) 


AfjycosO r 
(h+rr cofjd)* 


(h + rr cosd) r 


ie 


dp _ Sfjb. 


dr 


(H- 

t: 


4m rcosd 




(1+ y rcQsd 


g>r(i4:£S23®,» 


Integrating with respect to *r’ using tables of integration [33 


p _ 6 m h K 


I? r*cos*d 


1 + 


ycosdr 


2(1^. too?®!)' 


^ 4Mycogg 

•c® ® 

h r cos & 


In 


(1+ rco^j ^ 


1-H 


ycosdr 


2ii^ roo^^S f 


E 


1? 


£ ^ 2 + rcosdr/h ^a ^ ^ >y/E ^ 


1 + rcosdr/h 


1 + ycosdr/h 


+ (1.11) 


Where and are constants of integration. 
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1.6 Approximation 


Both equations (1.9) and (1.11) are much too cumbersome and an 
approximation, neglecting the radial variation oi film thickness, 
is considered and the pressure distributions so derived are 
compared with the results from eqs, (1.9) and (1.11). 

Thus film thickness is now given by. 


h = E + rr cos© (1.12) 

^ = >rr sine (1.13) 

ill 


and ^ = E + rx- cose ' (1.14) 

crt M 

Substituting in the Reynolds equation (1.1), with t = r 

in 

(h®^) = + IEm h + IZi^r^ji^cose 


3 


(h*^) = -6pyrjjjSiiie» + 12M(h + rr^j^cose) 


(1.15) 


1-7 Hydrodynamic Effect with Approximation 

From equation (1.15), the pressure distribution for only the 
hydrodynauJiic effect is given by 
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(1.16) 


“TP 


Intdffratlng with respect to ’r’ 


u» _ 


-6«rr^ sin® wr + C 

in 1 


Integrating again with respect to 'r’ 


-6Mrr 2 C r 

p = sir^ 


Where and are constants of integration 


i.e. p = - r_. sin^ w r* + ~ — + C 


C. r 


9 M 


With boundary conditions p = ps at r^ and p = 0 at r 


P- = 


- r_ sliO o.r' + ^ 


9 m 


+ C 

3 2 


and 0 = 


^ r. sine o>r^ 
1.3 m o 


C r 

+ -i— 2 + c 
1.8 2 


Substracting (1.18) from (1.19) 


-p = z M 

s ^3 ■m '”o "1' h* 


\j 

r sio? » (r* - r’J ) + — ^ (r - 
m 0 1,30 




te C = +7 r 


sindw (r* - r^) - p, 

l^s in o 1 s 


(1.17) 


(1.18) 

(1.19) 


(1.20) 
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and = 




■ 7 — = T ^ (rj + r ) 

(r^- r^) 


+ r sta9 w r* 

h* “ ° 


ie C, = - r sinS <■> r r, + p , ■ °-- - j 

* . » tt o 1 s tr r , ) 

h o i 


(1.21) 


Substitutinff for and in (1.17) 


^ r sti^ o> r* + -^ 3#irr^ sin©o> (r, + r^) 

IB ^ SI 1 O 


^ r sire « ps < °, ^ , 

h o i 


-3^^^ sia9 o> 


ie p, = r |r - (r +r. )r + r r . I 

*“ (E + rr_oos9)* 1- ° ^ ® 


^ (rt-t!) 


m 


■o "i^ 


1-8 Squeeze Effect with Approximation 

From equation (1.15), the pressure distribution for only the 
squeeze effect is £iven by 


ar 


/u* 

(h-;5j) 


IEm (h + r r_ cos®) 

III 


(1.23) 
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Intetfrattnff twice with respect to *r’ 


• C r 

= (E + rr„ coa0) r* + — ^ + C 

^3 In 


(1.24) 


where and are constants of integration. 

For boundary conditions, p = 0 at r = r^^ and p = 0 at r = r^ 
(The condition p = ps at r = r^ has already been incorporated 
the hydrodynajjiic pressure) 


0 = (h + r r^ COS0) r* + ■ + C 

h* “ ‘ h 


(1.26) 


and 


0 = 


C IT 

(h + r r CO50) r* + - + C 

,9 m « V* ^ 


(1.26) 


Subtracting (1.25) from (1.26) 


0 = -Si- (h + >- r. oosS) <r' 
j. m o 




i.e. C = -€iM (H + y r^ cosd) (r + r . ) 

3 ra v> X 


(1.27) 


in 
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and 


C = ^ {h + P r cost?) (r. r, ) 


rn 


o 1 


(1.28) 


Substituting for C and C in equation (1.24) 


6.w(h + r r cose* ) , 

m f 2 


sq 


h + j-- r 009^ )' 
ra 




(r- + r . ) r + r,r 
O 1 o 


J (1.29) 


1.9 Comparison of Results from the Complete and Approximate 
Equations for Pressure 

Numerical data used in the comparison are as below. 


^ = 0.1 Pa. s i.e. SAE 10 oil at 35‘C. 

r = 0.052m 

u 

r^ = 0.048m 

so that , r = 0.05m 

m 

h = 2.54 X lO'^m 


Since 


r r 
h 


• 0 . 1 , i . e . r 


max 


0.1 X 2.54 X 10 
0.05 


= 5.08 X 10 radians 


= 5.0 X 10 radians 


The shaft speed '&>’ is 3000 rpm 
i.e. CO = 50 rev/s. 


The value of r is also chosen as equal to 50 c/s. (This is just 
a convenient choice ; r does not have to be equal to co . ) 
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It may be noted that in each cycle the total traverse of r is 4 


times (and not 2rx ) . 

luolX 

r=4x(r ^ )x50 rad/sec 

rfi&x 

= 10*^ rad/s 

For tiie numerical comparison, we choose a tilt r, which is 
the maximum tilt 

ui-clX 


= 2.5 X lO"^ radians. 

Summary of Values 
(j - 0 . 1 Pa . s 

h = 2.5 X 10~^ ffi 

r - 2.5 X 10 radians 

r = 0.052 m 

o 

r^ = 0.048 m 

r = 0 . 05m 

m 

w = 3000 rpm ' 

r = 0.01 rad/sec 


half 
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eaa ]m mu mz ms <m 

ThetaCDegrees) 

riG. L5 HYDRODYNflHIC FILM PRCSSURE COMPflRISON 


I 


PressureCPaJ 
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1.10 Begultg 

Fig. 1.5 showg the hydrodynamic film pressure distributions, using 
the complete and approximate equations i.e. equations (1.9) and 

(1.22) respectively. Fig. 1.6, shows the squeeze film pressure 
distributions, using the complete and approximate equations ie 
equations (1.11) and (1.29) respectively. In both cases the error 
between the complete and approximate curves is below 2 % ( the 
numbers are given in Appendix B) and hence not visible on the 
graph. 

Thus comparison of the results from -the complete and 
approximated equations show negligible differences in p(d). One 
concludes that the variation of film thickness in the radial 
direction has a negligible effect on the overall circumferential 
film pressure distribution, and may be neglected. 


1.11 Combined Effect of Film Forceg 

Due to varying circiuftferential clearance in the seal interface, 
the fluid film has a converging and a diverging region. In the 
converging region, the viscous shear flow generated by the 
rotation forces the fluid into an ever - decireasing space, thus 
generating a positive pressure, which tends to force the seal rings 
apart against an external spring force. In the diverging region, 
a sub-atmospheric or negative pressure is developed. The negative 


27 





Fig. 1.7 Hy'drodynamtc film pressure distribution with and without 
large sealed pressure 
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FIG.- 1.8 FILH PRESSURE DISTRIBUTION 
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pressure, if present, leads to localised cavitation in the form of 
vaporisation of the fluid or release of dissolved «ases from it. 
The sealed pressure is assumed to he sufficiently high to ensure 
that the net pressures do not fall below atmospheric, thereby 
avoiding cavitation. 

Referring to figure 1.7, the constant part of the pressure, 
i.e. the sealed pressure p^, contributes only to the axial force 
on the ring, and not to the moment. If p is sufficiently large, p 
can be positive everywhere and only the sero-average (i.e. 
sinusoidal) part of the net p contributes to the moment. 

Figured. 8) shows the circumferential variation of 
hydrodynamic and squeeze film pressures. The hydrodynamic film 
pressure is anti-symmetric about the axis defined by d = 0* and 

180 and therefore contributes to a moment purely about that axis 
(see figure 1.9). The squeeze film pressure is symmetric about 
the same axis and therefore contributes to a moment purely about a 
perpendicular axis. Thus the moments generated by hydrodynamic 
and squeeze action are about mutually perpendicular axes. 
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Chaxvter 2 


2- 1 Dynamics of the Heehanlcal Face S«al 

A schematic figure of the mechanical face seal model beinff 
analysed is shown in figure 1.3. The rotor-seat is rigidly 
mounted on the shaft with no axial translating freedom and perfect 
alignment. The stationary seal ring is flexibly mounted on 
bellows. Therefore the stator ring has an axial translation and 
two rotations about mutually perpendicular diameters. Both faces 
of the seal are flat, and the basis for generation of film 
pressures is the angular misalignment (tilt) of the stator ring. 

The hydrodynamic and squeeze film pressure, and the spring 
forces, result in moments acting on the stator ring. The stator 
ring motion is analysed as a case of rigid body motion about its 
geometric centre. The geometric centre and mass centre of the 
stator ring are assumed to be coincident. Considering the tti.c 
rotational degrees of freedom of the stator ring about mutually 
perpendicular axes, the body has two equations of motion. 


2.2 The Co-ordinate system 

In the derivation of the expressions for pressure no 
refeirence was made to a bi-axial tilt of either ring. Instead, it 
was simply assumed that a net inclination y existed between the 
two rings and the film thickness was defined on that basis 
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[equation (1.12)}. If those results for pressure are to be 
directly usable, any subsequent choice of reference frame must 
preserve the conditions under which they were derived. Thus it is 
necessary to define a co-ordinate system in which the tilts can be 
described as about a single axis for all time. By definition, 
that axis must always be horizontal, i.e., parallel to the plane 
of the rotor ring. 

The co-ordinate system used in subsequent derivations is 
shown in figure 2.1. The X,Y,Z axes comprise an inertial 
reference frame, fixed in space with the origin at the geometric 
centre *0’ . The X-O-Y plane is parallel to the rotating seat and 
the Z axis is the centre line of the shaft. The x,y,z axes 
constitute a rotating reference frame. The x-O-y plane is always 
in the plane of the stator ring. The angle of tilt y, is the 
inclination of the axis of the stator ring relative to the 
vertical Z - axis and is also a measure of the angle between the 
planes containing the X,Y axes and the x,y axes. The angles r 
and ^ are the Euler angles which specify the position of the 
x-y-z rotating reference frame with respect to the X-Y-Z inertial 
reference frame. The x-y-z rotating reference frame precesses 
about the Z-axis with angular velocity y, such that the x - axis 
always remains in the horizontal plane (i.e. in the X-O-Y plane), 
and the y-axis always passes throuidi the point of maximum film 
thickness. Thus the x-y-z axes are the principal axes of the 
stator ring. The stator ring also nutates about the x-axis with 
angular velocity y . 
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Moaent. of tlo»en.tuM Kquationa 

The x-y-z axes and the stator ring have conunon nutation 
speeds about the x and y axes because the x-y plane is always in 
the plane of the stator. The x-y-z axes also have a spin about 
the z-axis which is not shared by the stator ring. 

Referring to fig 2.1, the absolute angular velocity w of the 
x-y-z axes in an inertial reference frame is seen to be 

« = Tjw + iw siny + ^ cosy (2.1) 

X y Sa 

The components of velocity are thus "x ~ ^ (2.2) 

Wy = V' siny (2.3) 

o = y cosy (2.4) 

z 

The absolute velocity of the stator is 
S = + “SSOg = iy + jysiny 


And the components of velocity are. 

Q = y 

(2.5) 


fty = y siny 

(2,6) 


= 0 

(2.7) 
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With principal axes of the present kind and different ansular 
velocities for the body and the rotating x-y-a frame, the Euler 
equations are [43 , 


“y "l * hz “t "y 


(2.8) 


n w + I Q 

y yy y zz z x xx 


w 
X y 


(2.9) 


M - ~ !», O &> + I O a> 

z sa a XX X y yy y x 


(2,10) 


Since the stator ring is axissmunetric, we shall use the ssrmbols 


I = I = I 
o XX yy 


( 2 . 11 ) 


I = I 


aa 


( 2 . 12 ) 


Substituting for the components of velocity from equations (2.2) 
to (2.4) and (2.5) to (2.7) in the Euler equations , 


"x ^ Ie ■ loV'^slny oosr 


(2.13) 


My = Iq (¥» ainy) + cosy 


(2.14) 


M„ = 0 

A 


(2.15) 


i.e. 


M„ = 


Ijr - siny cosy 
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(2.16) 


My. = 2I^ir COST + siny 
For small r. small tilts, 

sinr 55* r 

and cosy =» 1 (2.17) 

Using the relationships (2.17) in equations (2-16) 

{2.18) 

My, = 2I^vr + r (2.19) 

These are the Euler equations for the stator ring. 

The left hand sides of equations (2.18) & (2.19) represent 

moments acting on the stator ring. In this model, these moments, 
which are derived subsequently, are due to internal fluid film 
forces and the external spring force. 


2-4 Fluid Film Gfmerated tfcMsents 
The Film thicimess 

Bef erring to figure 2.2, the varying circumferential clearance in 
the seal interface (i.e. the film thickness) when viewed from the 
inertial reference frame is a right running moving wave. 
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Fis. 2.2 Model for film thi<^ess 
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Theiref ore , 


h = E+yr_ cos (^ -■ rt) 

n 

(2.20) 


(2.21) 

X* * * # # 

^ = n + yr^ij^ coa(^-lf't) + yr^^^ \^sin (^-yt) 

(2.22) 

pressure equations 



Substituting equations (2.21) and (2.22). in the Reynolds equation 

( 1 . 1 ) 


(b*^) = 6 m sin - vt) (2v^-«i) + 12 m (h 

+ cos(^ - yt) (2.23) 

la 

Using the same procedure as sections 1.7 and 1.8, the 
hydrodynamic and squeeze film pressure equations are given by. 


*Wd 


^r_sin(^-V't ) ( 2y^-» ) 

— ^ : 5- 

CE+yr cos(^ - ¥'t)] 

n 


(r - r) 

<V ■^i>' ^ ^ ®s (r^- r^) 

(2.24) 


P 


sq 


^{h + yr.cos{^-i't)) , 

— ^ i- Cr“ 

CEtrr^coaC^ “ vt)3 

HI 


< V Vi^ 


(2.25) 
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Fl«. 2.3 Model for celculattn* Bonexits 
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Tbo aoiBgt e^mattona 

Rdfonrins to fiffure (2.3), tho nomont about the X axis due to the 
hydrodynamic film pressure is given by 


‘“hy<i>X = J J r oo,* (r<*Klr)J 

r ^ 


(2.26) 


The limits of integration are r = r to r^ and ^ = 0 to 2rr . 
Substituting for from equation (2.24) and noting that the 

sealed pressure term is constant throughout the circumference and 
therefore does not contribute to any moments. 


^^d^X “ 3w'rjj^ (2r-«) 


I 


2- 2 
r Cr - 




rj^)r 


ViJdr 




cos^ 8in(^-yt)di^ 

Ch+ yr cos(-^-yt)* 
111 


(2.27) 


Referring to the appendix A and substituting the results (Al) 
(which defines a quantity R)and (A4) in equation (2.27) 


«Wd>X = 


”3|Myr^(2iy-») R w slort 

HI 


(2.26) 
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SinllftrlFf "th® ittOBi®n1; abou'b bh® Y axis, du® "to 'th® hsTdrodytianlc 
fluid film prassure ia. 


2n 
o 


‘“hjrd’v = J I r aiB» (rdpdrjj 


(2.29) 


Ti 0 


Subati-tuting for from aquation (2.24), (2.24) 


<“hy<l>y = r’cr‘-(vr^)r+ r^rj^ldr 


&r 

f 


aln^ ain (^-■yt)d^ 

CH + rg cos(^ - vt)]* 


(2.30) 


Raf erring to th® appendix A and subatitution the reaults (Al) and 
(A5) in e<iuatlon (2.30), we have 


^“hyd^Y 


3Myr_(2y-v>) R rr coay't 

m 


(2.31) 


Similarly, the moment about the X axis due to the squeeze film 
pressure is 


r 2nr 
o 

(Mg^)j^ “ J J (2.32) 

r o 
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The 8<iuee*;e film pressure contains a term in E. The quantity 'E is 
the time rate of change of the mean film thickness Ti and is 
non-sero when the seal closing force or external axial load varies 
with time. It may also be non-sero in the presence of transient 
film forces brought about by any perturbation in the 
geoBMtry. The face load determines the equilibrium values of the 
mean film -thickness, consistent with the amplitude of the film 
waviness. We shall assume here that the film profile represents 
an equilibrium configuration and that there is no externally 
induced vibration. Thus we may assume a constant H, for which 
E = 0. This assumption does weaken the model for the dynamics of 
the seal but is Justified on the ground that the present analysis 
is a first attempt at understanding face seal dynamics and the 
inclusion of the extra unknown Ei will introduce significant 
complexities in the model; for example, E cannot be constant as 
that would imply a monotonic departure from any postulated h. If 
E is included in the analysis. E must also be allowed for. 

Substituting for p ^ from equation (2.25) and neglecting the 
contribution of E to the squeeze - generated moments. 


2n 




= SMyr 


Wj Cr*-(V r^)r t r^r^lr'dr J ^ 


COS0 cos(^fr-yt)<l^ 
cos(^-^)3* 

nil 


(2.33) 
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Beferins the appendix A and substituting the results (Al) and 
{A6) in equation (2.33). 


(H 




Gnfjr^cog^t 

W 


(2.34) 


Sittilarly, the Boment about the Y axis due to squeeze film 
pressure is 




Srr^r^EJrsirm't 

11 


(2.35) 


The hydrodynamic moment which is purely about the y axis is given 
by the resultant of equations (2.28) and (2.31). .bus 


^“hydV " 


ftt^ir^R {¥'-w/2)r 

M 


(2.36) 


Similarly, the squeeze moment, which is purely about the x axis, 
is given by the resultant of equations (2.34) and (2.35). Thus 


(M 


sq^x 


enfjTjBy 

m 


(2.37) 
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i 


springs 


Fig. 2.4 Model for the Bellows 



Flg.2.5(b} Free body diagram of the stator ring 


showing spring forces 



2.5 Hc>iMm.1;g Itoe to the Bell.<wg 


The support arrazarenent oonslstlns of the bellows is modelled as 
equivalent of four equally spaced springs, each of stiffness 'k’ 
acting at the mean radius r^. The springs are assumed to lie on 
the X and y axes [fig. 2.4]. 

At any given tilt, which is always purely about the x axis, the 
spring force is distributed synoietrically about the y axis and 
therefore produces a moment purely about the x axis. Figure 
{2.5(a) shows the deflection of two springs on the positive y axis 
and negative y axis respectively, when the stator has tilted by an 
angle r about the x-axis. The positive x axis is in a direction 
out of the plane of the paper. 


Spring displacement is given by 

A - 1 - 

where, 

1 = 1 - A - h 

new 

6 = initial spring compression 

Substituting for h from equation (2.20), in (2.39) 

1^^„ = l-'S-H - yr. cos(^-vt) 
new m 


(2.38) 


(2.39) 


(2.40) 
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On the positive jr axis* ^ = vt, therefore 

1 „ = l-<5-E - rr^ (2.41) 

new n 

Substitutine equation (2.41) in (2.38), the spring displacement 
of spring *1’ is, 

A = A + E - rr^ (2.42) 

i. w. 

Similarly, on the negative y axis, ^ = ^t t n, and the spring 
displacement of spring ‘2’ is therefore, 

♦ 

A = d + E + rv^ (2.43) 

z m 

if 

Taking moments about the x axis, (fig. 2.5 (b) 

(*^sp) " -k(6+E+yr^)rj^ cost ^ k C‘5+E-rr||^)rj|^ cosy 

(**sp) ” "■2kr^coay (2.44) 

For small y, cosy Se 1, thus. 



-2kr* y 

w 



is equal to aero. 


(2.45) 
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2.6. Hfcft Iquattomi of Motloai 


Subatitutlna for the aiowsnte about the x and y axes given by 
eauatlona (2.37). (2,45) and (2.36) respectively, in the equations 
of »otion (2.18) and (2.19), the equations of motion are 

6rr#ir ^ * 

[giJ *■ 2krj|^ = (2.46) 


SfTAir R 

- - (¥'“<»>/2)y = 21 ^ + I V y 

U* o o 


(2.47) 


Let, 


A = 


- GftfJxJR 

Wk 

o 


Iji* 

o 


(2.48) 


and B = 2kr;^/I^ (2.49) 

thus, the equations of notion are, 
r -i- ir + (B-iy*)r = o 
and y + 4 Cl — ^ " ¥r/2i^)y = 0 

The equations (2.50) and (2.51), as they stand; have time 
dependent coefficients for r . In fact if y is a time dependent 


(2.50) 

(2.51) 
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quanti'ty, "the e^iuations are aon-linear and become intractable. 
The rotor however has a constant ansular velocity w. The 
precision ansular velocity of the rotating reference frame along 
the z axis, vis, w. occurs In response to « and may be assumed to 
be constant when o Is constant. Thus w Is assumed to be zero and 
equation (2.51) becoaes. 

f + i (1- = 0 (2.52) 

^ 2w 

In which w Is constant. 

Equation (2.50) is now a second-order, linear, ordinary 
differential equation in y. having constant coefficients. Its 
solution will be the solution for a damped spring-mass system; the 
damping term coefficient is positive because R is negative 
(eq. 2.48), and so the solution will correspond to the over 
damped, under-damped or critically damped case depending on the 
damping coefficient. Equation (2.52), on the other hand, is a 
first order, linear, ordinary differential equation which can have 

a tton-trivial solution of an exponential growth or decay. That 
is, equation (2.52) does not admit of an oscillatory solution. 

The quantity ^ represents a wobble of the stator ring. 
Although this quantity, has been assumed to be constant, it is an 
unknown and so it is appropriate to have the two equations (2.50) 
and (2.52) to describe the motion. It is necessary, however, for 
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these eaiiations to give the same solution for r tor all time. It 
follows that. If eqi. (2.52) does not have an oscillatory solution, 
neither can eq (2.50). A simultaneous solution to equations 
(2.50) and (2.52) is now sought. 


Differentiating equation (2.52) with respect to time, me get. 


jp + ^ (1- — ) y = 0 


(2.53) 


Cowpiirliig i5i<iuB.'fcloxi (2*50) wt'feli. (2.53) 

Ar + (B - ¥*)r = ■§ (1- -^)r 

^ 2V' 

i.e. -I (1+ -^) r + (B4®)y = 0 (2.54) 

^ 2r 

Equation (2.54) has an exponential solutions. Therefore let 


r = 


ro® 


xt 


(2.55) 


and 


r 


= xr 



(2.56) 


Substituting for r and r from equations (2.55) and (2.56) 
respectively, in equation (2.54) we get 


60 



(2.57) 


(B-y* I 
2r 

Substituting for i. from equations (2.57) in (2.55), the 
exponential solutions is 



r 



(B4^) 


1 + 


<•> 


2^^ 


(2.58) 


Three cases now present themselves, depending upon the relative 
siaes of and B. 

CASE 1. If V'* > B, r grows exponentially. 

CASE 2. If < B, r decays exponentially. 

CASE 3. If = B, y = i.e.y does not change. 


The 'CASE 3* is of principal interest. Using the definition of B 
(e.g. 2.49), it says that the face inclination r<«ialns constant 
if. 


w* = 2kr*/I^ 

m O 


i.e. V - ■/ -1^ 




A'iJ. 


(2.59) 

(2.60) 
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Equation (2.60) shows that the wobble frequency increases with 
increasing mounting stiffness and decreases with increasing 
transverse moment of inertia of the stator ring. This is as it 
should be. A "softer" support or a greater inertia will both 
make the statof ring more sluggish and its response frequency will 
be less. Conversely, if the stator ring has zero inertia, it 
will respond to the slightest impetus and will not sustain any 
moment . 

From equation (2.52), the corresponding condition for constant r 
is 



i.e. ¥' - "/2 (2.61) 

Substituting for ^ from equation (2.61) in equation (2.60) 



We now look for conditions which give other possible 
solutions to equations (2.50) and (2.52). 


Equation ( 2 . 50 ) also has an exponential solution. 
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Therefore using the same definition of r as in equation (2.55) the 
characteristic equation for (2.50) Is 

Kx + (B - V*) = 0 (2.63) 

where the subscript II refers to the second order equation (2.50) 


i.e. X 


/ 


- 4(B-y^) 


IX 


(2.64) 


Equations (2.64) indicates the conditions for real X^ (it must be 
real to simultfiuaeously satisfy eq. (2. 52)) is. 


a; 4 (B-V*) (2.65) 

Just as the ‘CASE 1' above indicate" that for y* > B , r grows 
exponentially, so also does equation (2.64) indicate a growing r 
for the same condition. 


Using the definitions for A and B from equations (2.48) and 
(2.49) respectively, 



( 2 . 66 ) 
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(2.67) 


Equation (2.67) sets a minimum value on y for real . 

Similarly, aquation (2.52) also has an exponential solution and 
the characteristic equation is. 



i.e \ = (l-«^/2y) (2.68) 

For both equation (2.50) and (2.52) to have a common solution, X^ 
and x^^ must be equal. Thus from equations (2.68) and (2.64). 

"I (l-w/2^) = ^ ± I / A*-4 (B-f*) (2.69) 

Simplifying and rearranging terms we get, 

+ {A_ ^ _ Q (2.70) 

The four roots of the quart ic equation (2.70) are, ( ppendix C) 
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and. 



B - 



+ 


1 

2 




Afo 

2 


(2.71) 


From equation (2.71), the condition for which v is real is. 


B > -4^ + .^ 

D - — T 2 


(2.72) 


i.e. B > I (| + w) 


(2,73) 


Using the definitions for A and B given by equations (2.48) 
and (2.49) respectively. 


2kr^ 6fT/ur„|Rj 

- Ifl ■ • 


m* 








+ c»> 


(2.74) 


Simplifying and rearranging terns, 


2k rl _ _ „ Sn^rj^jEj 






> 0 


(2.75) 


m 


Equation C2.75) sets a condition on 1?, namely that the minimum h 
is. 


1“ 


3rrAijR| 

4kr_ 

ni 


Cd 



(2.76) 
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Substituting for the half frequency solution i.e. equation (2.62) 
in equation (2.76), 


As H is a positive quantity, the only acceptable condition 
for li* is. 


« (1 


/ 7 ,] 


(2.77) 



3ti/j jRj 
4kr^ 


jw (1 + 



Using eq. (2.62) again for the half-frequency wobble 
condition, the above expression becomes 




1 . /T 


(2.78) 
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CHAPTER 3 


3.1. Main Results 

It is found that, witti a small angular misalignment between the 
seal faces, the radial variation in the film thickness can be 
neglected because it has almost no effect on the hydrodynamic and 
squeeze film pressure. Consequently, it also has negligible 
effect on the fluid film-generated moments and thus on th.e face 
seal dynamics. 

For the analysis of the face seal dynarciics, the mean film 
thickness is assumed to be constant. It is postulated that the 
stator ring wobbles in response to the rotation of the rotor and 
that, the shaft angular speed being nominally constant, the stator 
has a steady wobble frequency. The stator is flexibly mounted and 
is modelled as a flat but tilted ring. The wobble is modelled as 
a rotation of the stator tilt axis about the shaft centre line. 

With the above assumptions, two equations of motion are 
obtained for two mutually perpendicular components of the net 
tilt. It is found that tliese equations give only one meaningful 
solution, that of a steady-state wobble with constant amplitude 
as well as frequency. The corresponding wobble frequency is 
one-half of the shaft or rotating ring frequency. This finding 
would appear to corroborate the experimental observation of Burton 
and Etsion C53. Working with a set-up that duplicates t?ie present 
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Fi£. 3 



1 Siffnala showing a wobble frequency approximately 
half the frequency of rotation [5] 
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theoretical model, they found a steady operating condition in 
which the flexibly mounted stator ring wobbled at a frequency 
close to half the shaft speed. Their result is reproduced here in 
fig. (3.1), in which two oscilloscope traces are shown. The lower 
beam is a timing pulse, which se3rves to indicate the frequency of 
rotation of the rotor, whereas the upper beam indicates the wobble 
frequency. It may be readily seen that there are approximately 
half the number of pulses on the upper beam as on the lower beam. 

Looking back at the expressions for pressure, it is seen that 
the steady state wobble, for which the amplitude r is constant, 
leads to zero squeeze film pressure unless the mean film thickness 
varies with time (eq.1.29). This occurs because r - 0. 

Interestingly, the hydrodynamic pressure (eq. 1.22) also becomes 
zero because of the wobble frequency v' being one - half of the 
shaft speed except for the constant term in p^, the sealed 
fluid pressure. This constant term provides some axial force 
against face closure if p is significant, but does not generate a 
moment. So the steady wobble solution leads to zero moment. The 
fluid film still exists but has no resilience. It also generates 
no pressure with which to seal the leakage. In other words, the 
hydrostatic pressure gradient across the seal face, due to the 
sealed fluid pressure, causes fluid to leak out and the seal to 
fail. 
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The constant r or half -frequency wobble state is akin to a 
neutral state of equilibrium, with no fluid film resilience. Any 
disturbance of this state will therefore be unresisted by film 
forces. An increased or decreased tilt will be accompanied by a 
small departure in the threshold value of ^ (eq.2.60) and a shift 
to the exponential growth or decay solutions respectively. In the 
latter case, the stator tilt will be neutralised whereas in the 
former the stator and rotor rings will tend to come into physical 
contact. Either event will mean seal failure. 

The situation is analogous to the half -frequency whirl in 
hydrodynamic journal bearings. In that case, the journal centre 
whirls at one-half of the shaft speed resulting in zero 
hydrodynamic pressure and zero load-support in the bearing. The 
half ”f r'eQuency wobble of the seal face given by the present 
solution also results in zero fluid film pressure. 

The condition for h (eqn. 2.78) is a means of calculating 
the mean film thickness in order to estimate the leakage at the 
half frequency wobble condition. 


Other dynamic states are no doubt possible. To obtain them, 
a fuller analysis will be required. For example, if the 1 term is 
included in the analyses (i.e. if axial vibrations of the seal are 
allowed), that will give a non-zero squeeze - film moment. The 
extra equation required to solve for the extra unknown Tx(t) will 
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be the second order differential equation in Ti(t), describing the 
axial translatory dynamics of the seal ring. Based on the present 
experience, it can perhaps be anticipated that the said equation 
will also be non-linear. Of course, a complete dynamic analysis 
will have to accommodate a time-varying wobble ^(t) but that, it 
is seen, will lead to strongly non-linear equations of motion. 
Finally, the hydrostatic pressure gradient may also generate 
transverse moments in the flexibly - mounted seal ring [6] and 
influence its dynamic behaviour. 

The present half -frequency wobble solution requires two 
conditions to be satisfied simultaneously. These conditions are 
stated in eqs (2.60) and (2.62). Equation (2.62) is the necessary 
but not sufficient condition for eq. (2.61). In other words, the 
half frequency wobble does not necessarily follow if eq. (2.62) 
is satisfied. Thus it is suggested that the design of the seal 
should be done with an eye on conditions (2.62) and k and so 

chosen that, for the shaft speed dictated by the particular 
application, condition (2.62) is not fulfilled. Specifically, in 
order to avoid a case of exponential growth of y or a case of 
constant y, the condition. 



must be satisfied. Thus k and I must be so chosen for the shaft 

o 

speed dictated by the particular application. 
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CONCLUSION 


The flexibly mounted stator ring has two degrees of freedom. Its 
motion, in response to the rotation of the rotor ring consists of 
a wobble and is described by two non-linear differential equations 
for two mutually perpendicular components of the net tilt of the 
stator. The equations are reduced to linear form by assufliing 
constant wobble frequency and conditions are found for the tilt to 
increase, decrease or remain constant. It is found that when the 
condition for constant tilt is satisfied the wobble frequency of 
the stator ring is one -half of the shaft or rotating ring 
frequency. This state results in zero fluid film moments, causing 
failure of the seal. Equation (2.62) gives a relationship between 
the wobble frequency of the stator ring and the frequency of 
rotation of the rotor ring in terms of support stiffness, moment 
of inertia and mean radius of the stator ring. It is suggested 
that the design of the seal should be done with an eye on 
condition (2.62), and k and so chosen that, for the shaft speed 
dictated by the particular application, condition (2.62) is not 
fulfilled. Specifically, in order to avoid a case of exponential 
growth of r or a case of constant r , the condition. 


CO < 

o 



must be satisfied. 



APPENDIX A 


The lategral 
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O 

(r*- (r + r. )r*+ dr = f- -(r.+ r,)-^ + 

J o 1 0 1 I^D u 1 4 u X 3 J ^ 


R = -3" (!■’- r’)- (r^4- r. ) (r^- r*) t ^ r^r. (r*- r*) (Al) 

Where R is a negative constant. 

The Journal Bearing Integrals [7] 
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1 * 8ln(^-yt) ^ joo 

2«* _ Cl+« cos(^-Wt)3* * 



«he„, h 

(l-« ) ll+« oos(*-*^) 

SubstltutlM for in 



1 

e sin(^-^) 

« sin(^-^) -n ^ - il 

2e‘ 
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2 • 4 1 2 I 
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an 
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m 

o 

for small r> « 1. 
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Therefor© , 


Jtn 
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From (A2) and CA3) 
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for small r, y* « 1, ** « 1, therefore 
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APPENDIX B 


<l) Coaparison of coaplata and approx, valuaa of hydrodynaalc praasura 
<A) Coaplata 


ThataCDagraaa) 


. 0 

eo. 0 

40.0 

60.0 
BO . 0 

100.0 
1B0- 0 

140.0 

160.0 

160.0 
SOO . 0 

eso. 0 

240.0 

260.0 
EBO . 0 

300.0 

320.0 

340.0 

360.0 

CB) Approxlaata 


PraBsura(Pa) 


.0 

B 989.54 
17322.41 
24259 . 39 
2 B 94 B . 09 
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26166.53 
21603.05 
11921 .31 
. 0 

- 11921 .31 
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- 26186.53 
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. 0 
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.0 
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(fi) Comparison of coaplsts and approximats valuss of squsszs prsssurs. 
(A) Complsts 


Thstat Osgrsas ) f*rsssurs( Pa ) 


.0 

-66355.37 

eo.o 
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. 0 

£ 0.0 
ao. 0 
60.0 
80 . 0 
100.0 
1E0.0 

140 . 0 

160.0 
180.0 
£ 00.0 
££ 0.0 
£40.0 
£60.0 
£80.0 
300.0 
3£0. 0 

340.0 

360.0 


-6634£.7£ 
-6£681 .96 
-5£558.60 
-35658. £1 
-1E994.75 
13689.67 
41430. £3 
66145.35 
83377. £0 
89575 . 73 
83377. £0 
66145.35 
41430. £3 
13689.67 
-1 £994. 75 
-35658. £1 
-5£558.60 
-6£Bei .96 
-6634E.7E 


IZ 



AFFENBIZ C 


Solutions to the quart ic equation, 

¥* + (-^ - B) - (^)* = 0 
The quartic equation of reference £8], has the form, 

x* + ax* + bx* + cx t d = 0 

Comparing equations (Cl) and {C2), the coefficients are 

a = 0 

4 * 

b = - B 

c = 0 

d = - 

The resolvent cubic equation is of the form, 

7 * - hy* + (ac - 4d)y - a*d + 4bd - c* = 0 


(Cl) 


{C2) 


(C3) 


(C4) 


Substituting for the appropriate coefficients (C3), in equation 



/ - e-f- - B)r' + (^)'y - - B) (^)' = 0 




i.e ~ (~ ^~ 


(y* + (^)‘) = 0 


(C5) 


One root of the resolvent cubic eqn (04), required for solving the 
quartic equation is. 


y - T" 


(C6) 


Using the same nomenclature as given in the table being refered 
to. the required formulae are. 


R 


= / 


- b+y 


_ y 3a 


Zb t « 


= /IS 


o? ni. _ 4ab - 8c - a 
ET- 2b ^ 


(C7) 


Substituting the appropriate coefficients (C3), and one root of 
the iresolvent cubic equation (C6), in the set of formulae (C7), 
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R = y B 


_ A* . Afr) 

nr ^ ^ 


D = E 



(C83 


The formulae that give the roots of the quartic equation are. 


y 

A ,S 


- R + D 


y R + E 

• f* ~ 2 ? 


{C9) 


Substituting for D,E and R from equations (C8) 




(CIO) 
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